Effect of gauge boson mass on the phase structure of QED3 
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Abstract 

Dynamical chiral symmetry breaking (DCSB) in QED3 with finite gauge boson mass is studied 
in the framework of the rainbow approximation of Dyson-Schwinger equations. By adopting a 
simple gauge boson propagator ansatz at finite temperature, we first numerically solve the Dyson- 
Schwinger equation for the fermion self-energy to determine the chiral phase diagram of QED3 
with finite gauge boson mass at finite chemical potential and finite temperature, then we study the 
effect of the finite gauge mass on the phase diagram of QED3. It is found that the gauge boson 
mass m a suppresses the occurrence of DCSB. The area of the region in the chiral phase diagram 
corresponding to DCSB phase decreases as the gauge boson mass m a increases. In particular, chiral 
symmetry gets restored when m a is above a certain critical value. In this paper, we use DCSB to 
describe the antiferromagnetic order and use the gauge boson mass to describe the superconducting 
order. Our results give qualitatively a physical picture on the competition and coexistence between 
antiferromagnetic order and superconducting orders in high temperature cuprate superconductors. 
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Quantum electrodynamics in 2+1 dimensions (QED3) has been widely studied in the 
last years. It is well known that QED3 has two basic features: dynamical chiral symmetry 
breaking (DCSB) in the zero fermion mass limit and confinement l|-|8|. Moreover, QED 3 is 
superrenormalizable, so it does not suffer from the ultraviolet divergence which are present 
in QED 4 . Therefore QED 3 can be regarded as a useful tool by which one can develop 
insight into aspects of QCD. In parallel with a toy model of QCD, recently, as a low energy 
effective theory, QED3 is also of interest in 2D condensed matter systems, including high-T c 
superconductors |9l4l6l]. quantum Hall effect 17J and graphene 1814211]. 

At zero temperature, by solving a truncated system of Dyson-Schwinger equations, Ap- 
pelquist et al. first found that DCSB occurs only when the number of fermion flavors N 
is less than a critical number N c {3], 4]. At finite temperature, Aitchison et al. adopted an 
approximate treatment of Dyson-Schwinger equation for the the fermion self-energy in the 
1/N expansion 22| and found that the critical number of fermion flavors N c is temperature- 
dependent and chiral symmetry is restored above a critical temperature. Recently, Feng 



et al 



23- 



251 ] investigated the influence of finite chemical potential on the critical number 



of fermion flavors. It is interesting to study the influence of finite chemical potential and 



temperature on the critical number of fermion flavor simultaneous 
investigate the phase diagram of QED3 with gauge boson mass m a 



y. In this paper we shall 



261 ] at finite chemical po- 



tential and temperature. In particular, we study the effect of the gauge boson mass on phase 
structure of QED3. In addition, here we use DCSB to describe the antiferromagnetic (AF) 
order and use the gauge boson mass to describe the superconducting (SC) order (3, Q, 
and discuss the competition between m a and DCSB in chiral phase diagram of QED3 with 
gauge boson mass. We also note that the phase structure of QED 3 is also useful for the 
study of semimetal-insulator transition in graphene. 

In Euclidean space, the Lagrangian of QED 3 with N massless fermion flavors reads 



N 
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where £ contains the coupling between massless Dirac fermions and the U (1) gauge field, and 
£ is the gauge parameter. The gauge boson can acquire a finite mass m a via the Anderson- 
Higgs mechanism, which changes the original gauge boson propagator. The 4x1 spinor ipi 
represents the fermion field, the 4 x 4 7^ matrices obey the Clifford algebra {7^, 7^} = 26^, 
and i = 1, • ■ ■ , N are the flavor indices. For a real physical system, the number of fermion 



flavors N equals to 2 (the realistic model corresponds to N = 2, the two spin directions 



of the electron, see, for instance, 



28 



]). So, in this paper we investigate the chiral phase 
diagram of QED 3 with gauge boson mass m a at finite chemical potential and temperature 
for iV = 2 case. Adopting the rainbow approximation for the fermion-gauge vertex and 
neglecting the dependence of the gauge boson propagator on the chemical potential (this 
is a comm only used approximation in studying the dressed fermion propagator at finite \x 



23 



-25 



30M37j). we can obtain the Dyson-Schwinger equation for the fermion propagator at 



finite chemical potential \i: 



a 



d 3 k 



(2) 



where q = p — k, a = e 2 N, and D pu (q) is the gauge boson propagator in Landau gauge 



D,„{q) 



\V[IV 2 ) 



q 2 [l + U(q)}+m 2 a 

with 11(g) being the vacuum polarization function. 

At non-zero chemical potential, according to the treatments in Refs. 



(3) 



23 



24J , the inverse of 



the fermion propagator at finite /i can be obtained from the one at fi = by the substitution 
p^p = (p,p 3 + ifi): 



S 1 (p,fi) = S 1 (p) = ■ pA(p 2 ) + m(p 2 ), 



(4) 



where S 1 {p) = ry • pA(p 2 ) + m(p 2 ) is the inverse of the fermion propagator at \i — 0. The 
validity of the above formula (4) has been discussed in detail in Ref. 



Now, let us give a short review of some studies on the effect of the wave function renor- 
malization factor A(p 2 ). At zero temperature and chemical potential, the results in Eq. (2) 
show that when the 1/N order contribution to A(p 2 ) is included, the critical fermion flavor 
number takes almost the same values as the case where A(p 2 ) = 1 [4j. At finite tempera- 



ture temperature and zero chemical potential, a comparison of studies in Refs. [22|, |39| also 
suggests that the 1/N order contribution to A(p 2 ) only slightly changes the results qualita- 
tively. So, we expect that the the 1/N order contribution to A(p 2 ) is not important at finite 
temperature and chemical potential and we will take A(p 2 ) = 1 in this paper. Then, we can 
obtain the integral equation for the dynamically generated mass function: 



m(p 
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k 2 + m{k 2 ) 



(5) 



At non-zero temperature, for the purpose of obtaining a qualitative picture of dynami- 
cal picture of dynamical mass generation, we employ the following simplified gauge boson 
propagator ansatz 



40|, |41( in the presence of a finite gauge boson mass 

Q 2 + n (Q,/3) + m 2 



A^(g , Q, P) = ^ , J?*" , 2 , (6) 



where Ho(Q,(3) is defined as: 

n W,« = ^[^ + — exp(-— )|. (7) 

At finite temperature, the integration over the third component of a fermion loop momentum 
is replaced by an infinite sum over odd Matsubara frequencies. So the integral equation for 
the dynamically generated mass function at finite temperature and chemical potential is 
given by 

(R P) — V f rf2fe - m(P,ii,K) 

"HP. A*. ) NpJ^J (2ir)2Q2 + n (P,Q)+ml[K2+m%f3,ii,K) + (w n + i^y U 

where 

71 

P = (Po, Pj, ^ = |p|, Po = (2m + 1)-, 
k = {ko, k), A=|k|, fc = (2n + l)^, 

7T 

9 = (?o, q), Q = |q| = IP - k|, g = 2(m - n)-. (9) 
The sum over infinite Matsubara frequencies can be done analytically: 

OO 

J^(^ + ^) 2 + A 2 + m 2 (/3,/i,A) (10) 

1 ^ r 1 1 
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2 v /A 2 + m 2 (/3,/i,A) ^ n + ^A 2 + m 2 (/3,/i,A) ra„ — ^/A 2 + m 2 (/3 } fi, K) 
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2y / A 2 + m 2 (/3,/i,A) e^(/i - y'A 2 + m 2 (/3, /x, A)) + 1 e^(/i + ^A 2 + m 2 (/3, /i, AT)) + 1 
where we have made use of the formula 

y — - — = -^—- (ii) 

^u^-x e^ + 1 
Substituting Eq. (10) into Eq. (8), we can obtain the integral equation for the mass function 
with finite gauge boson mass at finite T and \i 

m[p,Li,P) = — - / a A . = 

8Nn2j ^K 2 + m 2 (/3, /i, A) x [Q 2 + n (/3, Q) + m 2 ] 

1 1 

(12) 



e^/i - J A 2 + m 2 (/3, //, A)) + 1 eP{fi + J K 2 + m 2 (/3, //, A)) + 1 



Here we also note that in graphene system the gap equation has a similar form 21] : 

1 1 1 , s 

x / = / = , 13 

e^(/i- v /K 2 + m 2 (/3,/i,K)) + l e/ 3 (/i+y / K 2 + m 2 (/3,/i,if)) + l 

where the Coulomb interaction V(Q,fx) = Tn &n f% QVjyi with the coupling strength A and 
the polarization function being given in Ref . 21] . In the context of graphene, the dominant 
interaction for Dirac fermions is the long-range Coulomb interaction, rather than gauge 
interaction. At finite chemical potential /i, however, the polarization function is proportional 
to fi as q tends to zero, which implies that the Coulomb interaction between Dirac fermions is 
statically screened and thus becomes short-ranged. The static screening factor \x appearing 
in the Coulomb interaction function in graphene system plays a similar role as the finite 
gauge boson mass m a in the effective QED3 theory of high-T c superconductor. For example, 
H suppresses semimetal-insulator transition in graphene system, while m a suppressed DCSB 
in QED 3 theory of high-T c superconductor. 
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FIG. 1. m(p 2 ) for several values of the gauge boson mass m a at N = 2, T = 0.0001 and \x = 0.0002. 

Now we employ the iteration algorithm to numerically solve Eq. (12). In the numerical 
calculation we set a = 1, thus every quantity with mass dimension is scaled and becomes 
dimensionless. In the real numerical calculation, we adopt an ultraviolet cutoff A = 10 5 , 
which is large enough to ensure that the calculated results are stable with respect to A. 
The calculated dynamically generated mass function m(p 2 ) is shown in Fig. 1 for N = 2, 
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T=1E-4 
H=2E-4 
m =0 




T = 0.0001 and // = 0.0002, and for several values of the gauge boson mass. From Fig. 
1 it can be seen that the dynamically generated mass m(p 2 ) is almost constant for small 
p 2 and tends to zero at large p 2 . This result is consistent with that given in Ref. [29]. In 
addition, we notice that m(p 2 ) decreases with the gauge boson mass increasing for any fixed 
value of p 2 . This reflects the following fact: Once the gauge boson acquires a finite mass via 
Anderson-Higgs mechanism , it cannot mediate a long-range interaction and suppresses the 
condensation of fermion-antifermion pairs. 

The phase diagram of QED 3 obtained in the presence of a small gauge boson mass for 
iV = 2 is plotted in Fig. 2. It is clearly seen that for each fixed value of the gauge boson 
mass the phase curve separates the DCSB phase (Nambu phase, where m(p 2 ) > 0) and the 
chiral symmetric phase (Wigner phase, where m(p 2 ) = 0). It can also be seen that the area 
of the region corresponding to the DCSB phase becomes small with m a increasing. This is 
understandable, because a finite gauge boson mass m a weakens the gauge interaction and 
suppresses the occurrence of DCSB. 
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FIG. 2. The phase diagram of QED3 for several gauge boson mass m a 

From Fig. 2, it can be seen that when m a is above a critical value m c ^' lt pa 0.018, the 
DCSB phase will rapidly disappear on the /i — T plane. This result gives a simple physical 
picture on the competition between the AF order and the SC order in high temperature 
cuprate superconductors, where the AF order dominates at low doping while the SC order 
dominates at higher doping. In addition, from Fig. 2 it is easy to find that DCSB disappears 

6 



and chiral symmetry gets restored when \x (or T) is above a certain critical value fi cnt (or 
T crit ), and the critical temperature decreases as the chemical potential \x increases. 

Just as was shown above, in QED 3 it is the dimensional coupling constant a (which is 
set to be 1 in our calculation) that sets the energy scale of the system. From our numerical 
results, it can be seen that the dynamically generated fermion mass, the critical temperature 
and chemical potential are quite small compared with this energy scale. Then, what is the 
reason for this? As far as we know, there is no simple way of seeing why the dynamically 
generated mass, the critical temperature and chemical potential are substantially less than 
the natural scale a. However, this is what seems to happen in all such calculations. For 

example, in the BCS calculation of the superconducting energy gap at T = 0, it is found 

i 

that the zero-temperature gap, A, is given by A « 2u£>e VN f , where lud is the Debye (cut- 
off) frequency, V is the volume and Np is the density of states. The quantity VNp plays 
the role of a dimensionless interaction parameter, and enters in a non-analytic way in the 
expression for A. This is essentially because the calculation is (strongly) non-perturbative. 
The natural scale is substantially modified by the non-perturbative exponential factor 
(for details, see, for instance, Ref. 42j). This exponential factor multiplying the natural 
scale is very similar to what happens in the T = QED 3 calculation. Something similar 
seems to be true at finite T as well. The fact that such calculations lead to dynamically 
generated symmetry breaking masses which are substantially smaller than a natural energy 
scale might help to explain the "hierarchy problem", namely the wide discrepancy between 
(say) Grand Unification scales and the weak scale (for instance, the idea that QED 3 might 
be a "laboratory" for studying the origin of different mass scales through non-perturbative 



dynamics was suggested in Ref. 



It is well known that the cuprate superconductor is a Mott insulator with long-range AF 
order at half filling. Once holes are doped into the Cu02 planes, SC order occurs and the 
long-range AF order disappears in the material. Recently, elaborate neutron scattering and 



scanning tunneling microscopy experiments 



44 



45] have found the competition between 



AF correlation and SC order in high-T c cuprate superconductor. It is clear that the gauge 
boson will become massive once SC appears in the cuprate superconductor. Here, in order 
to show the competition between the AF order and the SC order, following the approach of 
Ref. [271, we use the gauge boson mass to describe the SC order and use DCSB to describe 
the AF order. 
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To summarize, in this paper, in the framework of the rainbow approximation of the 
Dyson-Schwinger equation and employing the approximation of ignoring the \i dependence 
of the photon propagator and setting the wave function renormalization factor A{p 2 ) = 
1 , we obtain the gap equation of QED 3 at finite temperature and chemical potential. 
We investigate the phase diagram of QED 3 for several values of gauge boson mass m a 
at finite chemical potential and temperature for N = 2 case, which corresponds to real 
physical systems. In particular, we study the effect of a finite gauge boson mass on phase 
structure of QED 3 and find that a finite gauge boss mass m a suppresses the DCSB, which 
shows a competition between antiferromagnetic order and superconducting orders in high 
temperature superconductors. Besides, by a comparison between Eq. (14) and Eq. (15), we 
expect that the excitionic semimetal-insulator transition in graphene can be well described 
by the phase structure of QED 3 . 
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